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The average ground state energies for spin glasses on Bethe lattices of connectivities r = 3, . . . , 15 
are studied numerically for a Gaussian bond distribution. The Extremal Optimization heuristic 
is employed which provides high-quality approximations to ground states. The energies obtained 
from extrapolation to the thermodynamic limit smoothly approach the ground-state energy of the 
Sherrington-Kirkpatrick model for r — > oo. Consistently for all values of r in this study, finite-size 
corrections are found to decay approximately with ~ iV -4 / 5 . The possibility of ~ JV -2 / 3 corrections, 
found previously for Bethe lattices with a bimodal ± J bond distribution and also for the Sherrington- 
Kirkpatrick model, are constrained to the additional assumption of very specific higher-order terms. 
Instance-to-instance fluctuations in the ground state energy appear to be asymmetric up to the 
limit of the accuracy of our heuristic. The data analysis provides insights into the origin of trivial 
fluctuations when using continuous bonds and/or sparse networks. 



I. INTRODUCTION 

We study the ground state (T = 0) properties of spin 
glasses on Bethe lattices with Gaussian bond distribu- 
tion, which are also considered at low but finite temper- 
atures in Ref. [1]. In many ways, this study resembles 
that for the bimodal ±J bond distribution in Refs. [HQ. 
Yet, surprisingly, the behavior for finite-size corrections 
differ significantly between both distributions. 

Bethe lattices are r-regular graphs [J], i. e. randomly 
connected graphs consisting of N vertices, each having a 
fixed number, r, of neighbors [5j. We explore the large- TV 
regime of low-connectivity graphs, r — 3, . . . , 15, which 
are of great theoretical interest as finite-connected, mean- 
field models for low-dimensional lattice spin glasses [f| 0] ■ 
A great number of studies have focused on various aspects 
of this conceptually simple model to hone the complex 
mathematical techniques required to treat disordered sys- 
tems [El [sl— Tl3l| or optimization problems [T^4l9| . 

As before in Refs. [H, H, we use the Extremal Opti- 
mization (EO) heuristic |20h22| | to find approximations 
to spin glass ground states. In previous papers, we 
have demonstrated the capabilities of EO in determin- 
ing near-optimal solutions for spin glasses, the color- 
ing pro blem l2ll [23ll and the graph partitioning prob- 
lem [H, [H [2ji^ It is generally harder to find good 
approximations in complex energy landscapes with a lo- 
cal search heuristic, such as EO, for a problem with 
continuous weights (27| . but some encouraging results 
exist [28]. While our results appear to be sufficiently 
accurate for the prediction of energy averages for sys- 
tems up to size N — 2048, more detailed features, such 
as their fluctuations over the ensemble of instances (re- 
quiring higher moments of the energy), are less reliable 
for system sizes N > 256 and larger degree r. Hence, 
for Ising spin glass simulations with EO, discrete ±J 
bond weights are usually preferable. Based on the ex- 
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perience with the Sherrington-Kirkpatrick (SK) model, 
scaling properties of thermodynamic observables are gen- 
erally believed to be universal, independent of the de- 
tails of the bond distribution used. There, for instance, 
finite-size corrections are found to scale approximately 
with ~ 7V-2/3 [29M33|| (just as for the Bethe lattices with 
bimodal bonds) and even the average ground-state en- 
ergy density (esx) ~ —0.76317 ... is universal for any 
symmetric bond distribution [6j. It would be there- 
fore remarkable to find such distinct scaling behavior be- 
tween distributions at the level of finite-size corrections, 
as this investigation suggest. Some previous investiga- 
tions of Bethe lattices with Gaussian bonds [33l [34| have 
been consistent with JV -2 / 3 -corrections but were based 
on smaller sizes and significantly less statistics as in our 
study here. As a possible resolution, we would have to 
appeal to ad-hoc assumptions about higher-order correc- 
tions. (The distinct effects of discrete versus continu- 
ous bonds on defect energies in finite-dimensional lattice 
spin glasses have been studied numerically in Ref. (35| . 
and also with the renomalization group in Ref. [36[.) In 
turn, the variation of corrections with the details of the 
bond distribution may provide important clues towards 
extending replica theory to include finite-size effects. 

While the value of this work lies in exploring the range 
of finite-size scaling on sparse random graphs for spin 
glasses with differing distributions, it also provides a cau- 
tionary note about the analysis of data obtained from 
such systems. In a sparse system, such as a random 
graph or a randomly diluted lattice, (trivial) normal fluc- 
tuations in geometry [5l| or in the bond distribution can 
obscure the physical essence of the problem at hand. We 
will argue that in such a system we need to focus atten- 
tion to the actual "cost" C of disorder in terms of the 
frustration, instead of the energy E itself. Just consider 
a simple ferromagnet with fixed bonds J on an ordinary 
random graph Gn,p 01, or alternatively, with continu- 
ously distributed bonds J > on a fixed-degree Bethe 
lattice of size N . In either case, the ground state energy 
has all bonds satisfied, i. e. no cost in the number of 
frustrated bonds (C = 0), and ground state cost fluc- 
tuations exhibit a <5-peak, correspondingly. But by the 
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central limit theorem, the absolute sum of all bonds 

rN/2 

1=1 

has a normal distribution, inherited by the ground state 
energy fluctuations via E — 2C — B. We would claim 
that the relevant fluctuations at finite N are captured 
by C, not E, although the thermodynamic averages (E) 
and (C) are completely equivalent, of course. In any 
case, careful distinction is advisable in an environment 
of competing finite-size corrections, especially when av- 
eraging inherently finite samples in simulations. 

The effect of these trivial fluctuations is most pro- 
nounced in the study of fluctuations in the ground states 
of the systems, which have been intensely studied in re- 
cent years P, [29l - [33l [37l - l4l| . Our numerical data here 
shows that those fluctuations in the energy density would 
predict a seemingly interesting cross-over between non- 
trivial to trivial scaling between system size and degree. 
In contrast, the cost density fluctuations would predict 
a slow drift towards triviality at larger system sizes that 
is setting in earlier for larger system sizes and might be 
attributable to a decay in accuracy. As these results are 
somewhat inconclusive, in a related publication [42] we 
will show that these fluctuations on Bethe lattices with 
a bimodal distribution of bonds behave, again, similar to 
those of the SK model. 

In the following Sec. [Hi we introduce first the Bethe lat- 
tices we used in the numerical calculations. Then, we ad- 
dress the relation between energies and costs in Sec. IIIII 
In Sec. IIV1 we briefly describe the EO algorithm, which 
is amply discussed elsewhere [20-22]. In Sec. [V] we fi- 
nally present our numerical results. Some conclusions 
are presented in Sec. IVI1 

II. SPIN GLASSES ON BETHE LATTICES 

Disordered spin systems on random graphs have been 
investigated as mean-field models for low-dimensional 
spin glasses or optimization problems, since variables are 
long-range connected yet have a small number of neigh- 
bors. Particularly simple are Bethe lattices of fixed ver- 
tex degree r = k + 1 [Eli, which are locally tree- 
like with vertices imagined as possessing one up-direction 
and k downward branches. Yet, all r directions are fully 
equivalent in Bethe lattices, and there is no root vertex 
or any boundary. In comparison to the otherwise more 
familiar random graphs studied by Erdos and Reny [J], 
Bethe lattices at a given N and r avoid fluctuations in 
the vertex degree and in the total number of bonds. 

There are slight variations in the generation of Bethe 
lattices. For instance, to add a bond one could choose 
at random two vertices of connectivities < r to link until 
all vertices are r-connected. Instead, we have used the 
method described in Ref. [3] to generate these graphs. 
Here, all the terminals on the vertices form a list of rN 



independent variables. For each added bond, two avail- 
able terminals are chosen at random to be linked and 
removed from the list. Furthermore, for algorithmic con- 
venience, we reject graphs which possess self-loops, i. e. 
bonds that connect two terminals of the same vertex. 
Multiple bonds between any pair of vertices are allowed; 
otherwise it is too hard to generate feasible graphs for 
small N, especially at larger r. Since r remains finite for 
N — > oo, the energy and entropy per spin would only be 
effected to 0(1/ N) by differences between these choices. 

III. SAMPLING GROUND STATES ON A 
SPARSE GRAPH 

Once a graphical instance is generated, be it a Bethe 
lattice or any other sparse graph, we assign bonds Ji t j, 
here randomly chosen from a Gaussian distribution of 
zero mean and unit variance, to existing links between 
neighboring vertices i and j. Each vertex i is occupied by 
an Ising spin variable X{ G {— 1, +1}. The Hamiltonian 

H = - • / -. ( 2 ) 

{bonds} 

provides the energy of the system. For each instance I, 
the energy E^ is defined as the difference in the abso- 
lute weight of all violated bonds, C^ 1 ' (the "cost"), and 
satisfied bonds, i. e. 

E (D = C (D _ S (i). (3) 

the larger the satisfied bond-weight S^- 1 ' in the instance, 
the lower its energy E^\ While and vary de- 
pending on the spin configuration, for each instance 

£?« = CW + (4) 

is a constant, with B^ given by Eq. JT]). Hence, 

= 2C« - i?« (5) 

provides a direct relation between cost and energy of each 
instance. Thus, after averaging over a sample n\ of in- 
stances I of size N (denoted by (. ..)jv), we obtain by 
definition for the averages 

(E) N = 2(C) N -(B) N (6) 

and for the variances: 

a%(E) = W N {C) + a 2 N (B)-Acov N (C,B). (7) 

We observe that for the sparse-graph systems with con- 
tinuous bond weights under consideration here [similar 
to the ferromagnetic example in the Introduction, where 
a 2 N (B) = a% (E) ~ N and a 2 N (C) = cov N (C, B) = 0], 

°n(G) Z cov N {C 1 B)^a 2 N (B), (8) 

which we demonstrate in Fig. [T] Note that this condi- 
tion becomes less satisfied when graphs get denser, i. c. 
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r increases. That trend is indicated by arrows in Fig. [T] 
The situation described by Eq. © does not impact the 
relation between the averages in Eq. ((6]) in any way, in 
particular, not the thermodynamic quantities and their 
finite-size corrections discussed below. But it does sig- 
nificantly affect both, the error analysis for that data 
and the interpretation of the ground state energy fluctu- 
ations, each derived from the variance. Similarly a 2 N (B) 
can not be neglected in a fluctuating geometry, such as a 
random graph or a random-diluted lattice (43l . |44| , even 
if the individual bond- weights are sharp, \Jij\ = Jo with 
constant Jq. 

Eq. ([5]) does not apply for a discrete bond distribu- 
tion on a Bethe lattice [H, Q, where ajj (B) = 0; E and 
C are then equivalent stochastic variables in every re- 
spect. Furthermore, in a dense graph such as the SK 
model [45] , and independent of any symmetric bond dis- 
tribution, this situation is entirely distinct, as even for 
ground states almost exactly half of all bonds are vio- 
lated, i. e. 2C' / - ) ~ B^ for large N, and the energy 
appears merely as an extreme value within the normal 
fluctuations of 2C and B. In this case, using C to study 
fluctuations would be futile. As r is increasing, this trend 
is already visible in Fig. [1] the gap between the variance 
in E and C is closing, and is bound to cross over at some 
larger degree r. 

Eqs. ([7]|8]) imply that the standard error of the average 
energy is dominated by the error in B, 
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Figure 1: Plot of the variances and co- variances, appro- 
priately rescaled, in Eq. @ as a function of N for r — 
3,4,..., 10, 15. (Arrows indicate increasing r- values for the 
adjacent data sets.) All curves are either constant in N 
throughout, or appear to approach a constant. In this scaling, 
the variance for the bond weights B (red A) collapses for all 
r to a 2 (B)/(rN) = a 2 (| J|) /2 = § - ±, and the variances for 
the energies (green □) are of similar magnitude but declin- 
ing for larger r (top to bottom). At some fixed (larger) N, 
both, the covariance between B and C (blue o) as well as the 
variance of C (black o) alone, are quite small, but they are 
increasing with r (bottom to top). 
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which by Eq. ((HJ is larger than that of C for the range 
of degrees r used in this study. Therefore, our strategy 
for determining the average ground state energy in the 
thermodynamic limit will be based on an extrapolation 
for large N of the values for (C) N and the evaluation of 
Eq. © at N = oo using the exact value of the bond- 
density {b) aol here, 



of the form 



(c)n 



(C)n 
N 



(12) 



with a f=a (c)oo, even taking some higher-order correc- 
tions e(N) <C 1 into account. Additionally, we consider 
the fluctuations in the ground-state cost density, in par- 
ticular, the scaling of its deviation with finite size, 
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N~ 
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-= - 2 (\J[ 
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for Bethe lattices with a Gaussian bond distribution. 

More significantly, when Eq. ((HJ holds, the variance in 
the ground state energy fluctuations tracks the variance 
of B, making it apparently trivial (i. e. normal): 



'N 



(B) ~ N. 



(11) 



The scaling of <7jv (E) ~ A^~ p has been the focus of 
keen interest for various mean-field models recently P, 
[29l - [33l [37l - l4lT |. There, non-trivial behavior is typically 
associated with an exponent that obeys p > |.) We 
claim that due to Eq. (jlip any non-trivial deviation for 
sparse graphs would have to be found in the cost C, if it 
exists at all. 

In the following, we will therefore focus on the cost per 
spin c = C/N. We will study the finite-size corrections 



IV. r-EO ALGORITHM FOR BETHE LATTICES 

To obtain the numerical results in this paper, we used 
exactly the same implementation of t-EO to find ground 
states as in Ref. [2], except that we assign to each spin 
Xi in Eq. (HI) a "fitness" 
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<3> 



(14) 



Instead of using the local field of each variable, which may 
vary in range somewhat between variables due to the fluc- 
tuations in the bond values, we re-scale each variable into 
the same interval [—10, —9, —8, . . . , +10], thereby "hash- 
ing" the otherwise continuous state-space into a discrete 
set of up to 21 bins. Note that in this case the sum of 
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Figure 2: Extrapolation of the average ground state cost per spin for Bethe lattices of degree r = 3, 4, . . . , 10, and 15 on a 
l/7V"-scale for u = 1 (x), | (•), and | (□). The intercept of the fit with the ordinate represents the thermodynamic average 
(see Tab.lU, which varies weakly with the choice of U). For each r, a fit to the form a + b/N" , weighted by the error, consistently 
results in lo ~ | (see Tab. HTland Figs. 1314)) . Only rescaling to l/TVs virtually linearizes the extrapolation. 



the fitnesses is not proportional to the total energy, but 
good fitness sufficiently correlates with good costs for lo- 
cal search with EO to succeed, as explained in Ref. [20j . 

To evaluate the proposed r-EO algorithm, we have 
benchmarked over a number of exactly solved instances 
obtained with a branch-and-bound method. Such an ap- 
proach is clearly limited in attainable system sizes, here 
N < 64, and can only be executed for a small test-bed of 
instances due to the exponential computational cost of 
exact methods. For larger systems, we have also applied 
the r-EO algorithm to a small test-bed of 10 instances of 
size N = 2 10 for each value of r with 5-times more up- 
dates per run and compared results. In the worst case, 
about 30% of the instances for some r showed a system- 
atic error of about 0.1% or less. Since we sampled over 
ri] w 10 3 — 10 4 instances at this N, this systematic error 
is well below the statistical error of ~ 1/^/nJ > 1%. 

Alternatively, we can show that averaged proper- 



ties obtained with EO are in some sense self-consistent 
and/or consistent with certain theoretical predictions. 
For example, Fig. [S] shows the extrapolation for r — > oo 
of the already extrapolated thermodynamic (N — > oo) 
limit of the average energies for each r and N. Despite 
of its derivative nature, the EO data still reproduces the 
exactly-known ground state energy of SK very accurately. 

Finally, it should be kept in mind that settings which 
provide sufficiently accurate averages of a quantity may 
be less proficient in determining its higher moments, let 
alone its entire PDF. Ever higher moments are domi- 
nated by rare events (or the lack thereof for finite ni) 
ever deeper in the exponentially suppressed tails of the 
PDF. For instance, variances are based on (subtractions 
involving) higher moments and thus can be expected to 
have substantially boosted systematic errors compared 
to those quoted for averages. 
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Table I: Average cost per spin {c r ) N , plotted also in Fig. [2] The result for (c)oo is based on the value for a in Tab. HT1 



N 



Mi 



Mi 



Mi 



{C6)j 



Ml 



Mi 



Mi 



(cio)j 



16 
32 
64 
128 
256 
512 
1024 
2048 



0.09043(5) 
0.06934(3) 
0.05759(2) 
0.05087(1) 
0.04717(1) 
0.04491(1) 
0.04338(4) 
0.04302(2) 



0.20163(7) 
0.16825(5) 
0.14934(3) 
0.13842(2) 
0.13225(1) 
0.12853(1) 
0.12616(3) 
0.12577(4) 



0.3302(1) 

0.28476(6) 

0.25882(4) 

0.24387(2) 

0.23544(2) 

0.23026(2) 

0.22736(7) 

0.22666(7) 



0.4703(1) 

0.41254(7) 

0.37961(4) 

0.36103(3) 

0.35003(2) 

0.34363(2) 

0.33967(8) 

0.3384(1) 



0.6182(1) 

0.54844(8) 

0.50894(5) 

0.48600(3) 

0.47271(2) 

0.46465(3) 

0.4601(1) 



0.7732(2) 

0.6908(1) 

0.64416(6) 

0.61714(4) 

0.60153(3) 

0.59203(4) 

0.5863(2) 



0.9326(2) 

0.8383(1) 

0.78435(6) 

0.75324(4) 

0.73536(3) 

0.72440(4) 

0.7181(2) 



1.0976(2) 

0.9904(1) 

0.92878(7) 

0.89347(5) 

0.87325(3) 

0.86153(4) 

0.8545(2) 



1.9630(3) 

1.7954(2) 

1.6953(1) 

1.63804(6) 

1.60540(7) 

1.58710(9) 

1.5772(4) 



0.0420(1) 0.1236(1) 0.2236(1) 0.3352(3) 0.4542(5) 0.5798(8) 0.7103(9) 0.8459(9) 1.561(1) 



Table II: Fit of the data in Tab. U to (cr-)jv = a + 



r 


a 


b 




ndf x 2 / ndf Q 


3 


0.041995 


0.45 


0.81 


6 


92.5 







4 


0.123592 


0.70 


0.79 


6 


58.6 







5 


0.223676 


0.97 


0.80 


5 


63.6 







6 


0.335139 


1.22 


0.80 


5 


30.3 


6.3 x 10" 


-31 


7 


0.454128 


1.45 


0.79 


4 


38.2 


5.5 x 10" 


-32 


8 


0.579745 


1.72 


0.79 


4 


43.6 


1.3 x 10" 


-36 


9 


0.710283 


1.97 


0.79 


4 


40.4 


6.5 x 10" 


-34 


10 


0.845875 


2.32 


0.80 


4 


1.8 


0.12 




15 


1.560930 


3.65 


0.79 


4 


59.3 








Table V: Fit of the data in Tab. H] to (c r ) n = a+ -\ + 

N'S 



T 


a 


b 


c 


ndf 


X 2 /ndf 


Q 


3 


0.041941 


0.50 


-0.03 


6 


67.9 





4 


0.123668 


0.70 


0.01 


6 


62.9 





5 


0.223727 


0.97 


0.01 


5 


64.4 





6 


0.335286 


1.21 


0.02 


5 


33.7 


1.4e-34 


7 


0.454604 


1.37 


0.08 


4 


57.3 





8 


0.580229 


1.66 


0.07 


4 


59.0 





9 


0.710829 


1.90 


0.09 


4 


54.6 





10 


0.845843 


2.33 


-0.01 


4 


1.7 


0.14 


11 


1.561210 


3.44 


0.16 


4 


46.4 


4.4e-39 



Table III: Fit of the data in Tab. E to (c r ) N = a+^ + Tab i e V I: Fit of the data in Tab. Uto {c r ) N = a + jfe + § . 



r 


a 


b 


UJ 


c 


ndf x 2 / ndf 


Q 


r 


a 


b 




c 


ndf 


X 2 /ndf 


Q 


3 


0.041654 


0.61 


0.78 


-0.13 


5 


37.9 


5.5e-39 


3 


0.041500 


0.12 


0.62 


0.43 


5 


35.8 


9 x 10" 37 


4 


0.123219 


0.85 


0.77 


-0.12 


5 


50.1 





4 


0.123099 


0.32 


0.68 


0.49 


5 


50.5 





5 


0.223430 


1.06 


0.79 


-0.07 


4 


75.8 





5 


0.223414 


0.69 


0.75 


0.34 


4 


76.6 





6 


0.334402 


1.46 


0.77 


-0.20 


4 


20.2 


l.le-16 


6 


0.334213 


0.57 


0.69 


0.81 


4 


21.4 


1.2 x 10~ 17 


7 


0.452726 


1.86 


0.76 


-0.35 


3 


7.1 


8.6e-05 


7 


0.452095 


0.52 


0.63 


1.21 


3 


7.0 


9.9 x 10" 5 


8 


0.577866 


2.24 


0.76 


-0.45 


3 


7.6 


4.3e-05 


8 


0.576847 


0.54 


0.61 


1.54 


3 


5.6 


0.00073 


9 


0.708515 


2.47 


0.76 


-0.42 


3 


17.6 


2e-ll 


9 


0.707617 


0.73 


0.64 


1.61 


3 


15.0 


9.6 x 10" 10 


10 


0.845673 


2.38 


0.80 


-0.05 


3 


2.1 


0.1 


10 


0.845614 


2.02 


0.78 


0.35 


3 


2.0 


0.11 


15 




unstable 












15 




unstable 













Table IV: Fit of the data in Tab. Uto {c r } N = a+^- + ^-. Table VII: Fit of the data in Tab. Uto (c r ) N = a + -\ + 

N'S N 3 



r 


a 


b 


c 


ndf 


X 2 /ndf 


Q 


r 


a 


b 


c 


ndf xVndf Q 


3 


0.042773 


-0.21 


-0.45 


6 


2579.1 





3 


0.041656 


0.16 


0.37 


6 


34.3 


1.2 x 10" 


41 


4 


0.125082 


-0.37 


-0.74 


6 


1792.5 





4 


0.123011 


0.29 
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Figure 3: Plot of the \ 2 P er numbers of degrees of freedom 
(ndf) for a fit of the data in Tab. [Q to Eq. (O with only the 
parameters a and b over a range of fixed u), ignoring higher- 
order corrections (e = 0). Independently, for each r, the data 
shows a distinct minimum near 



V. NUMERICAL RESULTS 
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Figure 4: Plot of the fitted values for the finite-size scaling 
exponent u obtained in Tabs. [TT1 (□) and I VII (o) for the avail- 
able values of r. Based on an error of 0.01 estimated from the 
fits, all values from Tab.Hllare consistent with u> = | (dashed 
line), without trend for increasing degree r towards | (dash- 
dotted line). The values from Tab. IVll behave less uniformly 
but seem more centered around |. 



We have simulated Bethe lattices with the algorithm 
discussed in Sec. IIVI for r between 3 and 15, and graph 
sizes n = 2 l for I = 4, 5, 6, . . . , 11 to obtain results for 
ground state energies. Statistical errors in our averages 
have been kept small by generating a large number of 
instances for each N and r, typically nj ps 10 6 for N < 
256 and m 10 3 - 10 5 for N > 512. 



A. Average Ground-State Properties 

In Tab. HI we list the values of average costs per spin, 
(c) N . When plotted as a function of 1/N in Fig. [2j the 
average costs per spin for each given r clearly do not ex- 
trapolate linearly. Instead, we attempt a fit according to 
Eq. (|12|) with variable finite-size correction exponent co, 
at first ignoring higher-order corrections (e = 0). Listed 
in Tab. HU we find that for the whole range of connectiv- 
ities r studied here, the fitted scaling corrections appear 
to be consistent with u = 4/5, and a plot of the data 
in Fig. [5] on a rescaled abscissa, 1/iV 4 / 5 , produces linear 
scaling. In Fig. [31 we assess the quality of a fit, free of any 
assumptions about unknown higher-order corrections, by 
fixing the value of to over a range and then fitting for the 
remaining two parameters, a and b, in Eq. (|12p . Each of 
the independent data sets for r = 3, . . . , 10, 15 exhibits a 
strong preference for w « 4/5. In Fig. [H we plot those 
fitted values for uj, which within the range of degree val- 
ues r studied here suggest no trend toward the value of 
u>sk = § expected for the SK-limit r — > oo, irrespec- 
tive of bond distribution (3ll - l33l |38| . Plotting the data 
for w = | in Fig. [2] clearly does not produce a linear 
extrapolation as in Refs. 0, H[ for discrete bonds. 

In Tabs. IllltiVIH we have considered alternative fits to 



the data involving also higher-order corrections e(N) to 
the finite-size corrections, using the full form of Eq. (|12p. 
While there are some theoretical results (46|.|47|. obtained 
for the internal energy near T c in SK, possibly justifying 
u = | also below T c , little is know to higher order. The 
expansion for the free energy at T c in Refs. [!(| |47| pro- 
vided corrections of the form c\a(N)/N, which could be 
argued to eventually affect also the ground state energy. 
We have therefore attempted to fit the data also with 
higher-order corrections of that form, see Tab. Mil and 
simply c/N, another plausible form. [Unfortunately, a 
five-parameter fit with a higher-order correction of c/N a 
does not provide stable results for our data.] Assum- 
ing either correction provides the best-quality fit to the 
remaining four parameters, with the lowest x 2 value rela- 
tive to the remaining numbers of degree of freedom (ndf), 
shown in Tabs. Mil and IVII The former fit yields almost 
identical results to that without higher-order correction, 
Tab. mi only that the values for u> are consistently shifted 
down by a small amount, see Fig. @] The fit with 1/N 
drastically changes the fitted values. The values for ui are 
now somewhat consistent with |, but vary widely with 
degree r in Fig. U) Fixing lu = | improves the quality 
of fit (with one less parameter) for 1/N corrections to 
the best fit overall, see Tab. IVIII But it becomes en- 
tirely inconsistent with ]n(N)/N corrections, as shown 
in Tab. El unlike for fixed w = f in Tab. El 

In should be remarked that in all cases, the confidence- 
of-fit Q is essentially zero in light of the rather tight error 
bars obtained from the statistics for the average cost den- 
sities. We would argue that this is due to the inherent 
limitations of fitting this data down to small system size 
N with an asymptotic form valid for large N, Eq. (fT2"]) , as 



7 



-0.68 - 
-0.69 - 
-0.70 - 
-0.71 - 
-0.72 - 
-0.73 - 
-0.74 - 
-0.75 - 

-0.76,-.; 

-0.77 Ll 



□ Tab. II 

♦ Tab. VII 

linear fit, Tab. II 

linear fit, Tab. VII 



0.05 



0.1 



1/r 



0.15 



0.2 



0.25 



Figure 5: Extrapolated values of the energy densities (er)^ 
obtained from the cost densities (c r ) listed in Tab. U (yellow 
□). The data is plotted as a function of inverse degree, 1/r, 
and rescaled by a root of the degree, such that the large-r limit 
approaches the SK model. Similar to the corresponding plot 
m Ref. the data 

appears to reach the SK limit virtually 
on a linear trajectory: a linear fit (dashed line) to all data 
projects the ground-state of the SK model, (esx) = —0.76317 
(marked by x) as —0.76332, or to within an error of < 0.1%. 
An extrapolation of corresponding data from Tab. I VIII (blue 
o) misses the SK value noticeably (dash-dotted line). 



a stand-in for an entirely unknown function of N. Hence, 
these fits should not be dismissed on the basis of Q alone. 

The extrapolated values for a = (cr}^ obtained in 
Tab. [IT] are also listed in Tab. |TJ To demonstrate the 
quality of the extrapolation, we plot in Fig. [5] the derived 
values for the energy densities by way of Eq. (e^) = 

2 ( c r)oo - (Moo> usin S E q- (HDD- As in R-efs. @, [|, al- 
ready a linear fit to the extrapolated values reproduces 
the exactly known value of the SK model [|| [48|, |49j to 
within an error of < 0.1%, confirming to a high degree 
the numerical accuracy of the data. The obtained slope of 
the extrapolation, called /i in Ref. which is not ex- 
pected to be a universal quantity, evaluates to fi 0.36, 
much larger than for the corresponding problem with ± J 
bonds U 0|. The extrapolated costs in Tab. IVII1 also 
plotted as energies in Fig. [3J markably miss the SK value. 
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Figure 6: Scaling collapse for the deviation in the ground 
state fluctuations for the energy densities, crjv(e). The raw 
data for the Bethe lattices is plotted for the available de- 
grees r on top. On the bottom, the same data is plotted, 
now rescaled by the indicated powers of the degree, with the 
deviation also multiplied by \/~N such that trivial, normal 
fluctuations should plateau (as indicated by the continuous 
horizontal line). For fixed degree, systematic deviations ap- 
pear only for sufficiently small N < r 2 , consistent with a 
non-trivial scaling a(e) ~ N~ p at a value of p = | (dashed 
line) or even p = |, as expected for the SK model. If taken 
at face-value, one would conclude that for any fixed degree 
at large N merely trivial scaling is obtained and only the SK 
model (or any model with degree growing as r > v~N) has 
anomalous scaling, p > |. 



B. Ground-State Energy and Cost Fluctuations 

While the previous section has demonstrated the nu- 
merical accuracy of the data at the level of the average 
of observables, higher moments or a measure of the en- 
tire probability density function (PDF) of ground state 
energy fluctuations provide a more confusing picture. No 
clear conclusion can be reach on the basis of this data, 
even when heeding the implications of Sec. Mil In light 
of that discussion, we first present the (as we believe, 
incorrect) extrapolation for the ground-state energy fluc- 
tuations, followed by the corresponding discussion for the 
cost fluctuations. 



In Fig. we show the raw data for the deviations of 
the energy densities erjv (e) and their apparent collapse. 
All the data seems to approach trivial, normal fluctua- 
tions for sufficiently large system sizes, with a cross-over 
at ever higher degree r. Similar normal fluctuations for 
this model have been claimed by Ref. [33]. In fact, rescal- 
ing the data for a collapse indicates that the cross-over 
between system size and degree in the Bethe lattice oc- 
curs at N ~ r 2 . Taken at face value, this would be a re- 
markable result. While at any fixed, finite degree a trivial 
scaling is reached, only in the SK-limit, where system size 
and degree would scale direct proportionally, N ~ r, we 
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for small systematic errors based on the internal consis- 
tency of the average cost or energies, as is displayed, e. g., 
in Figs. [5] and [5l Though, this may prove insufficient to 
guarantee similar fidelity for higher cumulants, like the 
deviations in the cost. But if we assume sufficient ac- 
curacy for this data, it would imply that even for the 
cost deviations, like for the energies before in Fig. |51 ul- 
timately pure normal fluctuations may result either way 
when Gaussian bonds are considered. It would not be un- 
usual to find extended transient behavior in spin glasses 
with Gaussian bonds [501 • 



VI. CONCLUSION 



Figure 7: Plot for the deviation in the ground state fluctu- 
ations for the cost densities, o"jv(c). The data for the Bethe 
lattices is plotted for the available degrees, r = z, with the de- 
viation multiplied by yj N/r, such that trivial, normal fluctu- 
ations should plateau. No collapsible scaling regimes emerge, 
such as a non-trivial scaling <r(e) ~ N~ p with p = | (continu- 
ous line) and much less p = |, as expected for the SK model. 
Instead, all data appears to steadily approach trivial scaling 
at much larger system sizes N. 

would flow towards the left onto the non-trivial branch of 
the scaling in Fig. [5] Hence, the Bethe lattice results for 
Gaussian bonds (unlike for ± J bonds, as we will show in 
Ref. jUj]) are disconnected from the SK-limit. Although 
this may also explain the unusual finite-size scaling cor- 
rections to s» 0.8, as well disconnected from the SK-limit, 
we believe that this data collapse does not probe the 
true disorder-induced frustration. As we have shown in 
Sec. IIII1 energy fluctuations are dominated by the vari- 
ance in the Gaussian distribution of N bonds (although 
ever less so for larger r). 

Remarkably, when we plot the deviations <jn (c) for 
the cost density in Fig. a far more difficult-to-interpret 
picture results. Unlike for the energy densities in Fig. [HI 
there is no apparent cross-over but the data instead veers 
steadily towards normal fluctuations for much larger sys- 
tem sizes. While the data does not show any noticeable 
statistical errors, it is impossible to exclude a systematic 
bias in the sampling of ground states with the heuristic, 
that could manifest itself in a smooth drift away from 
any potential non-trivial scaling. We can only eliminate 
the systematic error up to N < 64 through comparisons 
with exact ground states obtained with a branch-and- 
bound algorithm. Above such sizes, we can only argue 



We have found surprising differences in the finite-size 
scaling behavior between a continuous, Gaussian bond 
distribution and previous results for a bimodal, ± J dis- 
tribution for spin glasses on Bethe lattices of degree 
r. While either distribution leads to equivalent results 
for thermodynamic averages that smoothly extrapolate 
to the exactly known SK results with identical scaling in 
r for r — > oo, the finite-size corrections term b/N u not 
only differ in the correction amplitude b but possibly in 
the scaling exponent u) itself. Only when higher-order 
corrections are postulated, more consistency can be ob- 
tained with the value w w | found for discrete bonds 
on Bethe lattices and for SK with either bond distri- 
bution. The value obtained here for Gaussian bonds, 
to « |, raises the question about an eventual cross-over 
to the SK value at higher r. No tendency toward such 
a cross-over is apparent in our study up to r — 15. In 
light of that, the vague expectation of some uncoopera- 
tive higher-order corrections to to = | seems preferable, 
but one might be forgiven to be struck by the solid per- 
sistence of u = | suggested by Figs. [5]and|U 

Our study of fluctuations in the ground-state proper- 
ties is not successful in determining clearly the scaling be- 
havior of the deviations. But it sends a cautionary note 
about the origin of fluctuations and the interpretation of 
data when simulating spin glasses on sparse graphs with 
continuous bonds or a randomly fluctuating geometry. 
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